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Abstract 



This article describes the investigation of morphological variations among two set of neuronal cells, 
Q ' namely a control group of wild type rat cells and a group of cells of a trangenic line. Special attention 

• is given to sigular points in the neuronal structure, namely the branching points and extremities of the 

dendritic processes. The characterization of the spatial distribution of such points is obtained by using 
a recently reported morphological technique based on forced percolation and window-size compensation, 
which is particularly suited to the analysis of scattered points presenting several coexisting densities. 
^ ■ Different dispersions were identified in our statistical analysis, suggesting that the transgenic line of 

neurons is characterized by a more pronounced morphological variation. A classification scheme based 
on a canonical discriminant function was also considered in order to identify the morphological differences. 

o . 

}Q • 1 Introduction 
O 

O ' During brain development, neurons form complex dendritic and axonal arbors that reach a characteristic 
' pattern and size PJ. The development of arbor shape is partly determined by genetic factors and partly 
O^' by interactions with the surrounding tissue (e.g. |21 El El)- Dendritic growth is known to be responsive 
J> , to various environmental signals, including synaptic activity and guidance molecules 013 El- The impact 
, of neurotrophins upon cortical neurons has gained special interest, since they are essential for neuronal 
^ I development, as well as for the maintenance of functional stability and plasticity in the adult nervous 
system 13 El Emm. 

Both neuronal activity and neurotrophins guide dendritic development via GTPase-dependent mech- 
anisms [51 ^2 ^1- Recent studies have drawn attention onto the small G-protein p21Ras in dendritic 
growth ^1 114L I15| . Like other GTPases, p21Ras regulates the phosphorylation of downstream kinases 
which trigger cascade mechanisms resulting in the regulation of enzymatic activities, ionic channels, cellular 
morphology and gene expression ^1 El ■ 

Effects of p21Ras have been shown to be essential for the normal functioning and plasticity of both 
the developing |19l 1201 1^ and the adult nervous systems (23 123 1^- Furthermore, enhanced expression 
of p21Ras is associated with neuronal restructuring after lesion or in the context of neurodegenerative 
diseases EHEiEIl- 
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investigated in p21H-Ras^'*'^^ transgenic mice |2HIEIEI- In this model the expression of transgenic p21H- 
Ras^*^'^^ starts postnatally around day 15, when neurons are postmitotic and the majority of synaptic 
contacts has been estabhshed. In par ticular cortical pyramidal neurons express the transgenic construct 
at high levels. The volume of the cerebral cortex of these transgenic mice is increased by approximately 
20% compared to wildtype mice. This increase has been attributed to the enlarged volume of the cortical 
pyramid al cells, but not to a raise in the number of neurons )29l l28[[TSlll4| . 

Dendritic morphology of pyramidal cells is very heterogeneous, not only in different areas (e.g. |3Up but 
also within the same region (e.g. jHlp. To set the population of neurons for investigation, only commissural 
neurons of layersII/III of the primary somatosensory cortex have been analysed in this study. By means of 
this restriction the analysis of a relatively homogeneous subpopulation of pyramidal neurons became possible. 




Figure 1: Retrogradely pyramidal neurons in layers II/III of the primary somatosensory cortex. Scale bar: 
50/iTO. 

In a previous study we found that enhanced p21Ras activity results in a dramatically enlarged dendritic 
tree. In both cortical layers II/III and V, the total surface area and the total volume of dendritic trees is 
greatly increased. This is mainly caused by increased dendritic diameter and by the appearance of additional 
segments [5^ . 

On the other hand, topological complexity of pyramidal neurons in layers II/III appeared hardly affected: 
ShoU analysis of both basal and apical dendrites revealed no differences between transgenic and wildtype 
mice regarding any parameters considered, i.e. numbers of intersections, branching points and tips |32|. 

These results suggest a rather proportional increase of dendritic tree size, without distinct changes in the 
space-filling properties. The aim of the present study was to substantiate these findings. 

In the present work we consider the use of a new procedure, known as percolation transform j!i8ll84| . which 
is particularly useful for characterization of spatial density of distributed points. This framework represents 
an enhanced alternative to the multiscale fractal analysis reported before pS]. The percolation transform 
procedure consists of two steps. The first is to induce an active percolation dynamics by performing an 
exact dilation of the points. During this process, the number of mergings between points and connected 
groups of points is recorded for each dilation radius. Having recorded the accumulative number of mergings 
as a function of scale, an essential normalization procedure is performed in order to compensate for sparser 
distributions which would have produced a much weaker signal otherwise. 

The overall enhanced sensitivity is a direct consequence of recording only mergings or shocks between 
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increase in area as considered in the (multiscale) fractal analysis This enhanced type of analysis 

is particularly suitable in the present case because the percolation transform approach is independent of 
size-related parameters like area and volume of the neuronal cells. We argue that changes in the global 
character of the percolation transform curves, derived from the reference points (i.e. dendrite extremities 
and branch points) of the dendrites of pyramidal neurons due to transgenic activation of p21Ras in the 
primary somatosensory cortex of mice, correlate with changes in the complexity of neuronal morphology. 

2 Materials and methods 
2.1 Biological 

We used data derived from experiments with three transgenic mice aged nine months (see ^ EHl for the 
establishment of the transgenic mouse line), as well as with three wildtype mice of the same age. All 
experimental procedures on animals were carried out in accordance with the European Council Directive of 
24 November 1986 (86/609/EEC) and had been approved by the local authorities. All efforts were made to 
minimise the number of animals used and their suffering. 

Mice were anaesthetised with Hypnomidate® (Janssen-Cilag) (lml/40g body weight intraperitoneally) 
and positioned in a stereotaxic apparatus. Prior to incision 0.1ml Xylonest® (Astra) was injected below 
the skin, the left parietal bone was partially removed. Using a Hamilton syringe for unilateral injection, 
Biotinylated Dextran Amine (11 20% BDA, 10.000 nominal molecular weight(MW), Molecular Probes) was 
mechanically delivered into the corpus callosum 1mm lateral to the midline and 0.5 mm caudal to the 
bregma, 1mm deep from pial surface according to the stereotaxic atlas of Frankin and Paxinos (1997). 
Seven days after surgery the animals were intracardially perfused in deep anaesthesia, first with saline (0.9% 
NaCl) for 1-2 min and then with a fixative containing 4% paraformaldehyde in 0.1 M phosphate buffer (pH 
7.4) for 30 min. The brains were postfixed in the same solution overnight, immersed in 30% sucrose for 
another 24 hours, and sectioned (160m) in the coronal plane on a cryostat. The free-floating sections were 
extensively washed in 0.05 M Tris Buffered Saline (TBS, pH 7.4), and then reacted with the avidin-biotin- 
peroxidase complex (1:100 in TBS, VECTASTAIN® Elite ABC kit. Vector Laboratories, Inc.) at room 
temperature for two hours. BDA was visualised by using 3,3-diamino-benzidine (Sigma, 0.025%) intensified 
with nickel-ammonium sulphate (Merck, 0.05%) in the presence of 0.001% hydrogen peroxide, diluted in 
TBS. Sections were mounted onto gelatine coated glass slides, dehydrated and covered with DPX (Fluka, 
Neu-Ulm, Germany). 

The retrogradely labelled pyramidal cells were reconstructed using Neurolucida® (MicroBrightField, 
Inc.), see Figure ^ The system allowed accurate tracing of the cell processes in all three dimensions and 
continuous adjustment of the dendritic diameter with a circular cursor. A motorized stage with position 
encoders enabled the navigation through the section in the xyz axes and the accurate acquisition of the 
spatial coordinates of the measured structure. All visible dendrites were traced without marking eventual 
truncation of smaller dendritic sections. This may have led to certain underestimation of the dendritic 
tree, especially in transgenic mice with a larger dendritic arbour. To gain both optimal transparency for 
optimal tracing facilites and at the same time a possibly complete neuronal reconstruction, sections of 160 m 
thickness were used. Thicker sections allowed only ambigous tracing of thinner dendritic branches. Shrinkage 
correction (300%) was carried out in the z axis, but not in the xy plane, because shrinkage was negligible in 
these dimensions ( 10%). 

The morphology files created by Neurolucida® need to be edited and converted to a simpler file format 
before they can be used for further analyses. We used for this purpose the freely available program Cvapp 
Cvapp is a cell viewing, editing and format converting program for morphology files, see Figure |21 It can 
be also used to prepare structures digitized with Neurolucida® software for modeling with simulators like 
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to the SWC format describing the structure of a neuron in the simplest possible way |[{7| . 





Figure 2: Pyramidal cells rendered with Cvapp. Displayed are random examples of each transgenic neuron 
(ceh SE8, left) and wildtype neuron (cell WT19, right). 

Each line encodes the properties of a single neuronal compartment. The format of a line in a SWC file is as 
follows: n T X y z R P. In turn, these numbers mean: (1) an integer label (normally increasing by one from one 
line to the next) that identifies the compartment, (2) an integer that represents the type of neuronal segment 
(0-undefined, 1-soma, 2-axon, 3-dendrite, 4-apical dendrite, etc.) , (3)-(5) xyz coordinate of compartment , 
(6) radius of compartment, (7) parent compartment (defined as -1 for the initial compartment). The sample 
of neurons edited, converted and analyzed (see below) comprised 28 of wildtype neurons and 26 of transgenic 
neurons. 

2.2 Computational 

The concept of the percolation transform [23 El is described in the following. As mentioned in Section ^ 
this method uses of a mechanism known as Minkowski dilation. Given a set of points 5' in i?'^, its Minkowski 
dilation by a radius r corresponds to the union of closed balls of radius placed over each point in S, and the 
volume of the dilated set is represented as V{ri). In order to ensure enhanced precision and computational 
efficiency, only those radiuses corresponding to the viable distances in the orthogonal lattice, the so-called 
exact distances |3H], are taken into account. 

Consider now that the set of points S is generated from a Poisson distribution with density 7 inside 
a cubic window of side L. Let C{di) be a function that counts the number of mergings between growing 
clusters at a specific distance di = 2ri as the whole set of points undergoes a dilation by radius . In order to 
allow scale uniformity we make the change of variables ki — hi{di). A graphic representation of this function 
is shown in Figure O for two different values of density 71 and 72. One can see that there is a characteristic 
scale for each density, where most of the merging dynamics takes place before reaching a plateau. The total 
number of points in S is given by jiV or 72 y as illustrated by Figure |3 The idealized scenario where this 
dynamics occurs is represented by the unit step functions. In our model, the position of such a transition is 
taken as the characteristic scale for the set S. 

Next, we need to relate the caracteristic scales to the densities of the generated Poisson set of points. We 
argue that the distance where all points suddenly colapse into a unique cluster is in fact the mean nearest 
neighbor distance. Knowing the nearest neighbor distance distribution |39| for the three-dimensional space, 
the mean nearest neighbor distance as a function of the density of the cloud can be expressed as 



0.554 



(1) 
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The smooth curve represents the behaviour of the shock dinamics for a Poisson model, while the unormalized 
Dirac deltas stands for the idealized model. Note that the rate of change of the cumulative number of mergings 
for higher scales tends to decrease fastly. As we intend to use these characteristic scales to detect multiple 
densities coexisting in the same set of points, we need a normalization procedure in order to enhance the 
signal for sparser distributions of points, therefore assigning the same importance to all scales. To implement 
such a normalization we take as reference the maximum value of the rate of change for one specific reference 
scale kj. The function we need should always give the same value irrespectively of the characteristic scale 
we are dealing with, as in the following equation 

Q(fc,) = :^e(3S-3^^). (2) 

This normalization procedure is the last step in the percolation approach. In the case of an experimental 
cumulative function obtained from a Poisson simulation with specific density. Equation El will produce the 
same intensity for any density value. When dealing with an experimental cumulative function involving 
an unknow mixture of densities. Equation |2] will yield detect the multiple characteristic distances actually 
occuring in the sampled space. 

We illustrate the potential of this procedure regarding coexisting densities of spatially distributed points 
for the following experimental example. Figure [S] presents an example pattern (inset) and the corresponding 
percolation transform curve. The example pattern consists of two Poisson clouds with distinct densities. The 
percolation curve for this example shows clearly the two peaks corresponding to the characteristics scales 
associated with the two different densities. Note that the difference in height of those peaks reflects the 
relative area occupied by each clowd. 

In the present work we consider the spatial distribution of the singular points of the dendritic trees, i.e. 
dendritic branch points and extremities. Such an approach (see j4ni I41| ) has shown to be an effective route 
when dealing with subtle geometrical changes while considering a single morphological/physiological class of 
neuronal cells. 




Figure 3: The smooth curve shown corresponds (schematically) to the accumulative mergings as a function 
of the characteristic scale for a Poisson distribution of points. The unit-step functions shown come from 
idealized assumptions. 



5 



Figure 4: The derivative of the cumulative merging function, as a function of the characteristic scale. 
Schematically as before the smooth curves stands for Poisson clouds and the Kroncckcr functions, represented 
as upward arrows, are associated with the idealized model. 

3 Results 

Figure IHl shows the resulting percolation transform curve for the two selected cells of Figure |21 The overall 
profile of these curves were the same, showing two preeminent peaks at clearly distinguishable scales ^. 
Although we can see that these two single cells presents a clearly distinc global features in their signature 
curve, this property is masked by the statistical variabitities of the whole set of cells, as shown by the scatter 
plot in Figure[3 This scatter plot is an example of various quantile-quantile graphs we produced using global 
properties of all percolation curves, such as those exemplified in Figure El For this purpose, we select a few 
global properties, namely the mean, the variance, the maximum value, the scale at which the maximum value 
occur and the monotonicity index [1^. The scatterplot of Figure is defined by the monotonicity index and 
the maximun scale, as an example of the overall variability of shapes. Taken pairwise none of these global 
measures leads to clear separation of the two types of cell involved in this study. 

Before ruling out the possible existence of a geometrical characteristic that may be correlated with the 
genetic treatment, we performed a more thorough statitical analysis of the data produced by the percolation 
approach. These consist of a principal component analysis and a cononical discriminant analysis j43j . Fig- 
urelHlshows the density profile for the principal component variable, showing a similar mean but a pronouced 
diference in dispersion of the considered geometrical features of the genetically treated cells. Tabled shows 
the performance of the those selected global features in providing a discriminant function for the two type of 
cells. While there is visible potential for classifyng one type of cells, namely the SE group, the discriminant 
function has a marginal performance concernig the other group, misclassifying almost half of them. This 
can be graphically visualized in Figure |51 





SE 


WT 


Error 


Posterior. Error 


SE 


20 


7 


0.2592593 


0.3891051 


WT 


9 


12 


0.4285714 


0.4767072 


Overall 






0.6878307 


0.4274310 



Table 1 : The result of classical discriminant analysis for the measures considered in the scatterplot of Figure[7| 
From this cross-validation table we observe a somewhat better classifacation rate for type SE while for WT 
the estimation is indecisive. 



Note that the scale axis is the logarithm of the actual distance among the points. 




Figure 5: An example showing the resolving ability of the proposed procedure for the pattern show at the 
inset. The two coexisting densities are clearly detected by the peaks occuring at their characteristic scales. 



4 Discussion 

The present study is one of the few which uses 3D data on neuronal morphology to quantitatively char- 
acterize the complexity and scaling properties of different neuron types. The morgphology of pyramical 
cells was investigated in two sets of neurons, i.e. wildtype and p21H-rasVall2 transgenic mice. The results 
obtained using principal component analysis indicates that the transgenic neurons have a greater dispersion, 
as revealed by the density profiles shown in Figure |S| This finding may suggest that the enhanced p21Ras 
activity in transgenic mice may lead to greater variety of the cell morphological phenotype. 

These results must be valuated with care, however. We know from other studies (e.g . Fernandez 
et al. 1994) that a measure like fractal dimension for example is not sufhcient to differentiate between 
cell classes. Alone it does not completely specify a cells morphology but increases classification accuracy 
as an additional parameter for morphological classification among several other parameters typically used 
for placing nerve cells in different classes such as soma diameter, maximum dendritic diameter, number of 
branches etc. The same situation may have occured here with the percolation transform approach. Although 
our findings suggest that transgenic mice pyramidal neurons exhibt a more variable dendritic morphology 
than the corresponding wild type neurons, this problem must be investigated in more detail. Clearly such 
study will include a more extended statistical analyses of the neuronal morphometrical parameters. 
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Figure 6: Example of the percolation curve for the two types of neuronal cells presented in Figure |2| 



g 0.6 
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Maximum scale 



Figure 7: A scatter plot defined by the monotonicity index (vertical axis) and the scale (horizontal axis) at 
which the maximum of the percolation transform curve ocurs. 
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Figure 8; The density profile of the principal component score for each class in this experiment. One can 
see a very similar mean and broader dispersion in the WT cells shown by the dotted line. 
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Figure 9: The graphical result of canonical analysis, showing the amount of discrimimation possible for both 
classes of neurons with the one dimensional variable. 
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